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1.  NTRODUCTION 

Gravity  waves  in  water  surface  are  characterized  by  a  small  steepness  value.  It  makes  possible  to  apply  the  weak- 
turbulence  approach  in  order  to  find  out  the  problem  solution  with  the  help  of  kinetic  equation  for  energy  spectmm. 
This  equation  was  formulated  by  K.Hasselmann  (1962,1963)  and  V.Zakharov  (1968).  In  terms  of  wave  action 
spectrum  the  equation  is  as  follows: 

+ki  -k2  -A:j)6(co+cq  -co^-a:^)x 

x[N ^{N  +  N ^)-  N^N{N^  +  N^)]dk^dk^dk^  (1) 

where  =N{k-  )  is  a  spectral  density  of  wave  action;  T (kjc j,k2,k^)  is  a  kernel  function  of  non-linear  wave 

interaction;  6  (k)  and  S(G))  are  the  delta -functions  describing  a  resonance  interaction  between  four  wave 
components: 


k  +A:^  =  k2  +k2  ;  (2a) 

CO+COi  =C02+W3  (2b) 


For  gravity  waves  —  ,  the  value  T(k,kj,k2,k^)~k^  is  a  homogenous  function  of  the  third  order.  An 

explicit  expression  for  T{k,kj,k2,k^)  can  be  found  out  in  original  papers  (Hasselmann  1962,  1963;  Webb  1978; 
Zakharov  1999). 

The  equation  (1)  has  been  a  subject  of  numerical  modeling  for  almost  three  decades  (Webb  1978,  K.Hasselmann 
1985;  Resio  and  Perrie  1991;  Polnikov  1993;  Komen  et  al.,  1994;  Komatsu  and  Masuda  1996,  Zakharov  1999, 
Lavrenov  1998,  2001  etc.).  However,  some  basic  properties  of  this  equation  are  not  clarified  till  now. 

Usually  the  equation  (1)  is  considered  to  preserve  standard  constants  of  motion  -  wave  action,  energy  and 
momentum: 
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(3) 


N  =  ^N{k)dk 


E  =  \()i{k)N(k)dk 
K=^kN{k)dk 


Actually,  only  is  a  real  constant  of  motion.  The  energy  E  and  the  momentum  K  are  only  “formal”  integral  of 
motion,  “leaking”  in  the  area  of  very  large  wave  numbers.  The  high  frequency  truncation  of  the  Hasselmann  equation 
leads  to  a  leakage  of  the  wave  energy  and  momentum  to  a  high  frequency  range,  whereas  the  wave  action  flux  is 
mainly  directed  to  a  low  frequency  range.  It  can  be  presented  as  follows: 


dN  ^  dE 

- =0;  - 

dt  dt 


dK 

dt 


=  -M; 


(4) 


where  P,  M  are  fluxes  of  energy  and  momentum  in  k  -space  directed  to  high  frequency  area. 

Preservation  of  energy  and  momentum  fluxes  to  high  wave  numbers  is  due  to  formation  of  weak  turbulent 
Kolmogorov  spectra.  The  main  theoretical  point  of  the  kinetic  equation  (1)  is  a  description  of  the  stationary  equation 
solution  of : 


(5) 


The  simplest  Kolmogorov  weakturbulent  stationary  solution  of  the  equation  (1)  was  obtained  in  1966  (Zakharov  and 
Filonenko,  1966): 


5((0)=ao 


^4/3pl/3 


(6) 


where  S[(X>}  is  the  energy  spectrumdefmed  by  the  relation: 


S((0)d(a  d^  =(0(k)N(k)  dk  =  — ^A^(co,  p  )r/co  d^ 

8 

In  general  terms  the  Kolmogorov  spectra  are  anisotropic: 


(V) 


5(co^|3)= 


4/3  pl/3 


CO 


gM 

coP 


(8) 


where  F  =  F(^,  |3)  is  a  function  of  two  variables. 

A  solution  (9)  averaged  over  the  angles  is  very  close  to  energy  frequency  spectra,  obtained  in  many  laboratory  and 
field  experiments  (Toba  1972;  Donelan  et  al.  1985). 

Now,  a  question  arises,  how  can  t  he  exact  Kolmogorov  solution  of  the  stationary  equation  (5)  approximate  solutions 
of  the  nonstationary  equation  (1),  especially  in  cases  of  applying  forcing  and  damping  sources?  It  should  be  noted 
that  for  the  first  time  an  attempt  to  give  an  answer  to  this  question  is  undertaken  in  paper  (Komen,  Hasselmann  and 
Hasselmann,  1984).  They  show  that  even  in  the  presence  of  input  and  dissipation  there  is  the  spectrum  solution 

5  ~  l/co"'  .  The  nonlinear  flux  is  found  out  to  be  generally  rather  strong,  so  that  relatively  small  deviations  from  the 

solution  S  ~  l/co"^  are  sufficient  to  generate  divergent  fluxes,  which  can  balance  nonzero  input  and  dissipation 
source  function  in  the  cascade  region  (Komen  et  ah,  1994).  Unfortunately,  at  that  time  it  was  difficult  to  get  an  exact 
estimation  of  basic  parameters  of  problem  solution..  A  stationary  solution  obtained  by  solving  nonstationary  equation 


(1)  should  be  developed  with  reliable  accuracy.  That  is  why  it  is  necessary  to  use  sufficiently  accurate  and  fast 
numerical  algorithm,  not  available  at  that  time. 

It  should  be  noted  that  a  full-scale  experiment  on  numerical  simulation  of  the  equation  (1)  is  performed  recently 
(Pushkarev  et  al.,  2002).  The  approach,  based  on  another  numerical  algorithm  (Webb  1978;  Resio&Perrie,  1991), 
existing  for  more  than  two  decades,  is  used.  As  this  problem  seems  to  be  rather  complicated  for  numerical 
simulations  and  the  results  are  obtained  for  the  first  time,  there  appears  a  necessity  to  produce  independent 
estimations  in  order  to  verify  the  above-mentioned  results  with  the  help  of  the  another  algorithm. 

This  is  a  motivation  for  writing  present  paper.  Numerical  results  of  the  equation  (1)  solution  taking  into  account 
wave  energy  forcing  and  damping  are  obtained  in  this  paper  with  the  help  of  the  numerical  algorithm  elaborated 
recently  by  one  of  the  authors  (Lavrenov  1998,  2001).  The  algorithm  is  based  on  numerical  methods  of  highest 
precision.  It  is  of  high  accuracy  and  calculation  speed. 


2.  PROBLEM  FORMULATION 

The  equation  of  spectral  wave  action  evolution  can  be  written  as: 

d  N 


dt 


■yiv  =G„,  +F 


(9) 


G..  ic 


where  the  non-linear  energy  transfer  function  (1);  Y  is  a  damping  increment  depending  on  the  frequency 

CO  as  follows  : 


Y  = 


Cl  w sin '  ((1  -  co/co„i„ )  (tc  /2)),  if  (a  < 
0,  tJ(o„i„<(o<co^ 

C2Wp(w/w^-l)',  ifcoxo^ 


(10), 


where  Cj  2  are  positive  constants. 

The  value  F  is  an  external  active  force:  F  =  f  N  .  The  function  /  =  /(CO,  P)  is  a  value  not  equal  to  zero 
within  the  frequency  range:  <CO<COf  .  It  is  equal  to  the  following  angular  function: 


|2cos"(Ap),iJ  cos(Ap)  >0  and  co„i„  <co<co^ 

f  =  \  ,  (11) 

[0,  in  another  cases 

where  2  is  a  normalizing  function,  providing  the  same  integral  value  for  the  various  parameters  n  and  A: 


% 

|/(P,(0)dp  =  C3(0y  (^2) 

-71 

where  C3  is  a  constant. 

So,  a  problem  is  posed  in  such  a  way  that  the  whole  frequency  range  is  divided  into  sub -ranges  with 
different  energy  sources.  The  low-frequency  damping  domain  is  located  in  the  range  (CO  <  )  in  order  to 

stabilize  fluxes  directed  into  the  low  frequency  range.  The  energy  pumping  domain  is  located  within  the  range 
( <  CO<  03^  ).  One  of  the  most  interesting  frequency  ranges  is  the  domain  (00^  <  CO  <  CO^  )  without  any 
damping  or  pumping,  i.e.  it  is  the  so-called  “transparency  window”,  where  the  spectmm  is  formed  only  by  non-linear 


energy  transfer.  The  Kolmogorov  constants  should  be  defined  using  numerical  results  in  this  range.  The  high 
frequency  damping  is  located  in  the  domain  ( 03  >  03^ ).  Values  of  the  frequencies  ;  COy  ,C0^  are  defined  in 

such  a  way  that  the  appropriate  frequency  ranges  include  sufficient  number  of  grid  points  to  obtain  reliable  integral 
estimations  defined  by  spectmm  values  within  the  corresponding  frequency  ranges.  The  following  values  are  used: 

o^in-0.5;  cOy  -1.0,(Op-6.5 

In  order  to  comply  with  the  requirements  of  the  kinetic  equation  (1)  application,  the  coefficient  Q  in  (12)  is  defined 
as  ^3  “  0.001  satisfying  the  conditions  of  smallness  of  the  growth  rate  with  respect  to  the  corresponding 

frequency:  f  «CO.  The  coefficients  C12  as  well  as  the  appropriate  frequencies  CO^;„;COy,CO^  are  defined 

experimentally.  The  conditions  of  effectiveness  of  fluxes  absorption  and  minimization  of  the  pumping  and  damping 
intervals  with  respect  to  appropriate  frequency  are  used.  It  should  be  noted  that  exact  values  of  the  constants  Cj  2 

(satisfying  the  conditions  Cj  2  <  1  -0  )  are  not  so  principal  to  obtain  a  qualitative  solution  within  the  range 

«lnin  <  CO<  COp  .  The  values  Q  2  are  defined  in  such  a  way  that  the  energy,  wave  action  and  momentum  fluxes 

are  absorbed  by  dissipation  within  appropriate  frequency  ranges  and  by  numerical  accuracy  of  fluxes  estimations. 
Numerical  results  are  estimated  in  96  directions  and  50  frequencies.  Such  a  detailed  angular  resolution  is  used  for 
obtaining  an  accurate  estimation  for  narrow  angular  distributions  of  external  force  approximation  (12)  considered  in 
paper.  In  this  case  the  angular  increment  is  equal  to  A(3/27r  ~  0.065,  being  the  same  order  of  frequency  increment: 

A(0/(0~  0.068  .  Such  discritization  is  optimal  for  numerical  integration  in  general  case,  when  solution  is  not 
known  before  hand  (Lavrenov,  1998). 

As  the  solution  is  to  be  obtained  for  a  large  time  scale  (up  to  100000  seconds),  an  optimal  numerical  algorithm  of 
non-linear  energy  transfer  computation  is  used  (Lavrenov  1998,  2001). 

3.  NUMERICAL  RESULTS 
3.1  Spectrum  evolution 

Now,  the  numerical  results  for  the  isotropic  case  with  n  =  0.0  and  A=1.0  should  be  considered.  The  frequency 
spectrum  can  be  defined  as  a  function  of  wave  action  as  follows: 


S(C0)  = 


71  n 

|s(C0,p)cfp=  I  W  (CO,  p)2C0‘‘/g  "  dp 

-71  -71 


(13) 


The  spectmm  (13)  for  the  following  time  steps:  100,  300,  500,  1000,  3000,  5000,  10000,  30000  and  50000  sec., 
respectively,  is  shown  in  the  logarithmic  scale  (Fig  1). 

Two  different  stages  can  be  defined  in  the  wave  spectrum  evolution.  On  the  first  stage  the  spectrum  growth  is 
observed  within  the  range  of  the  external  force  impact:  0.5  <03  <1.0.  The  spectmm  is  quickly  increased  in  more 
than  6  orders.  Duration  of  this  time  interval  is  estimated  as:  t  =  10000  sec.  The  spectrum  becomes  almost  stationary 
at  f  >  30000  sec.  Its  values  coincide  exactly  with  the  spectrum  at  f  ~  50000  sec.  The  initial  and  final  stages  of 
spectral  evolution  are  presented  in  Fig.  la.  Within  the  range  0.8  <  CO  <  8.0  ,  the  stationary  spectrum  is  very  close  to 
the  Zakharov- Filonenko  one: 


^^'v“4.0±0. 1 

5  (CO)  =CO 


(14) 


Spectra  evolution  within  the  intermediate  time  interval  1 1000<  t  <  15000  sec.  is  presented  in  Fig. lb. 


a) 

- 1 -  t=100s 

- 0 -  t=300s 

- B -  t=500s 

- @ -  t=1000s 

S  {m^ s  /  rad)  ^ -  t=3ooos 

- ^ -  t=5000s 


Frequency 


b) 

- 1 - t=1 1 0OOs 

- $ -  t=1 2000s 

2  ,  B -  t=13000s 

S{m  si  rad)  0 -  t=i4ooos 

- A -  t=1 5000s 


Frequency 


a;  (rad/s ) 


CO,  (rad/s  ) 


Fig.  1 

Frequency  spectrum  evolution  within  time  intervals: 

a)  100<r<  50000  sec  ;  b)l  1000<  t  <  15000  sec  ;  ( — )  -  approximation  S  ~  C0^° 


3.2  Estimation  of  spectral  front  propagation 


In  wave  spectmm  the  evolution  of  two  effects  can  be  observed:  a  relatively  slow  growth  of  spectra  in  the  forcing 
range  and  a  fast  “front  propagation”  to  area  of  high  wave  numbers.  It  is  a  self-accelerating  process.  In  paper  (Komen 
et  ah,  1994),  this  concept  is  used  to  explain  high  frequency  spectmm  formation  and  the  migration  of  the  spectral  peak 
towards  lower  frequencies. 

A  formation  of  the  stationary  spectra  occures  within  a  finite  time  interval.  At  the  same  time  the  wave  energy 
propagation  to  high  frequency  range  is  observed.  The  speed  of  this  propagation  can  be  estimated  as  a  shift  of  the 
spectral  frequency  front  into  high  wave  numbers  (Fig.  2).  The  spectral  frequency  front  is  defined  as  a  high  frequency 
(“cut-off  frequency”),  with  spectral  density  value  being  one  order  smaller  in  comparison  with  the  limited  spectrum 

solution  for  a  given  frequency  within  the  range  (0^  <  (0  <  (0^  .  The  same  effect  is  described  in  paper  (Donelan  et 

ah,  1985). 


Fig.2 

Propagation  of  spectral  frequency  front  ^front  into  high  frequency  range 
1  -  data  approximation  (16),  2  -  numerical  data; 

According  to  theoretical  estimations  ^akharov  et.  al,  1992),  the  front  propagation  is  approximated  using  a  self¬ 
similar  solution  of  the  equation  (1): 

5(CO,f)  =  (fo  (15) 


Hence,  the  “cut-off  frequency”  is  increased  explosively: 


(16) 


^  front  — 


COo 


Numerical  estimations  obtained  in  our  experiment  are  in  agreement  with  these  theoretical  results  (Fig. 2).  A  self¬ 
similar  front  propagating  solution  can  be  explained  as  a  result  of  energy  diffusion  in  frequency  space  in  order  to 
reach  a  system  stable  state. 


4.  ESTIMATION  OF  KOLMOGOROV  S  CONSTANTS 


4. 1  Definition  of  Komlogorov’ s  constants 

According  to  the  weak  turbulence  theory  the  general  Kolmogorov  spectrum  is  defined  by  fluxes  of  the  wave  energy 
P  ,  wave  action  Q  and  momentum  M  (Zakharov  et  ah,  1992).  The  spectrum  is  assumed  to  be  symmetrical  with 
respect  to  the  reflection  P  — >  —  P  .  In  a  general  case,  the  spectrum  is  as  follows: 


5(0i|3): 


4/3pl/3 


CO 


cog  gM 
P  ’coP 


,COS  P 


(17) 


In  this  paper  only  a  case  with  no  wave  action  flux  at  infinity:  2  —  0  is  considered.  A  similar  solution  was  obtained 
analytically  (Zakharov  and  Zaslavskii,  1982)  for  the  infinite  frequency  interval,  when  energy  pumping  was  located  in 
the  low  frequency  range  and  energy  damping  in  the  high  frequency  one. 

For  the  large  value  CO,  the  function  F(f,  p  )  can  be  expanded  into  the  Taylor  series  for  small  ^  =  g  M  /coP  : 


P(^,  P)  ~  ttp  +  OC,  COS  P 


(18) 


For  spectrum  S (CO,  P)  ,  it  can  be  presented  as  follows: 

g  M  COS(  p  )  ^ 


5(C0,p)  = 


CO 


ttg  -i-aj  ■ 


V 


CO  P 


(19) 


J 


where  OCq  and  CXj  are  the  first  and  the  second  Kolmogorov  constants,  which  are  coefficients  of  the  spectral  density 
expansion: 


OCq 


1 


CO 


1  } 


71 

j5(C0,p)r/p; 


‘  K 


1 5(co,P)cos(P)  r/p 


(20) 


(21) 


The  energy  flux  P  ,  directed  into  the  high  frequency  range,  is  estimated  as: 

n  oo 

P=  J  jY(CO,P)5(C0,p)r/C0r/p 

-It  COp 

The  momentum  flux  is  estimated  similarly: 

71  oo 

M  =  J  Jy((0,P)5((0,P)— cosPii(oc/p 

-71  COp 


(22) 


(23) 


4.2.  Isotropic  case 


The  value  CXg  can  be  determined  in  isotropic  case  as  follows: 


_  5'(co)co'^ 

g4/3pl/3 


(24) 


According  to  the  weak  turbulence  theory,  the  value  OCg  is  a  constant,  not  depending  on  CO.  As  our  estimates  show: 
CXg  =  0.303+ 0.033  for  time  moment  t>20000  sec. 


4.3.  Non-isotropic  case 

Similar  numerical  simulations  of  the  non-linear  spectmm  evolution  are  fulfilled  for  non-isotropic  source  function 
generating  wave  energy  (12)  with  different  values  of  the  parametersn  and  A. 

The  time  evolution  of  wave  integral  parameters  (i.e.  total  energy,  wave  action  and  momentum)  is  shown  in  Fig. 3.  A 
full  stabilization  of  these  parameters  is  observed  at  f  ^  10000  sec. 


-e- 


Fig.3 

Time  evolution  of :  1-  total  energy,  2-  wave  action;  3  -  momentum 


Estimation  results  of  the  limited  spectmm  at  t=20000  sec.  for  the  most  narrow  angular  distribution  are  presented  in 
Fig.4  (frequency  spectmm)  and  Fig. 5  (frequency-angular  spectmm).  It  should  be  pointed  out  that  the  main  features 
of  the  frequency  spectmm  are  nearly  the  same  as  in  the  isotropic  case. 


S  {irP' s !  rad)  - 1 -  i 

- 2 


Frequency  (rad/s) 


CO,  (rad/s) 


Fig.  4 

Limited  frequency  spectrum  in  non -isotropic  case  ,  1  -  numerical  data;  2-  approximation  S  ~  (£) 
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Fig.5 

Limited  frequency-angular  spectrum 


As  for  non-isotropic  cases  there  is  an  additional  possibility  to  estimate  not  only  the  first  Kolmogorov  constant,  but 
also  the  second  one  using  (21).  The  corresponding  results  are  presented  in  Table  1 

An  angular  spectmm  width  D  (Fig. 7)  can  be  determined  as  a  proportion  of  the  frequency-angular  spectrum  value  at 
the  main  direction  to  the  frequency  spectrum  as  follows: 


D(a))=5((0,|3j/5(co) 


(25) 


D 


(0,  rad/s 


Frequency  (rad/s) 

Fig.  6 

Function  D((0) 

1-  numerical  results;  2  -  local  approximation:  D  ~  0.235  • 

The  function  D  is  shown  in  Fig.  6.  Its  maximum  value  is  achieved  at  the  frequencies  (0=  0.5  -F  0.6  rad/s.  After 

that  it  is  quickly  reduced  as  soon  as  the  frequency  is  increased,  showing  that  the  function  of  angular  energy 
distribution  becomes  wider. 

Final  estimation  results  of  the  Kolmogorov  constants  for  different  values  of  external  angular  functions  are  presented 
in  Table  1.  As  it  is  seen  the  estimations  of  the  first  and  the  second  constants  are  varied  weakly  depending  on  the 
external  angular  function  distribution. 


ESTIMATION  OF  KOLMOGOROV  CONSTANTS  FOR  DIFFERENT  ANGULAR  DISTRIBUTION  FUNCT  IONS 


TABLE 1 


Angular  function  distribution 
of  energy  generating  force 

tto 

Oi 

Isotropic 

0.303  ±  0.033 

COS  (P/2) 

0.308  ±  0.020 

0.218  ±  0.015 

cos^(P) 

0.324  ±  0.021 

0.239  ±  0.023 

cos"'  (p) 

0.31 1  ±  0.023 

0.242  ±  0.034 

It  should  be  noted  that  the  results  presented  in  Table  1  are  in  general  agreement  with  (Pushkarev  et  ah,  2002),  in 
which  the  following  estimations  are  obtained:  for  isotropic  case:  0.35  <  OCg  <0.45  and  for  non-isotropic  one: 
0.33  <ao<  0.37  and  0.18 <0(,  <0.27. 

The  discrepancies  between  the  above  mentioned  results  and  ours  can  be  explained  as  follows.  In  (Pushkarev  et  ah, 
2002)  there  are  some  differences  between  the  values  OCq  for  isotropic  and  non-isotropic  cases,  whereas  the  values 

OCq  should  be  constant  according  to  its  definition.  It  should  be  noted  that  these  estimations  are  obtained  within  a 
more  wider  confidence  interval  and  only  for  two  angular  distributions.  In  our  computation  numerical  simulations  for 
four  angular  distributions  are  produced:  from  isotropic  to  a  very  narrow  one  typical  for  wind  wave  generation. 
Another  reason  of  the  discrepancies  can  be  explained  by  the  fact  that  there  is  no  evidence  that  fully  stabilized 
numerical  solution  is  obtained.  (Pushkarev  et  al.,  2001).  Moreover,  our  estimations  can  be  considered  to  be  more 
accurate  due  to  the  fact  that  a  more  precise  numerical  algorithm  and  more  fine  angular  resolution  are  used. 

5.  COMPARISON  WITH  TOBA  SPECTRUM 


Toba  (1973)  made  careful  measurements  of  the  spectrum  tale  and  found  out  that  the  following  approximation 
describes  a  spectrum: 

sm  =  8gu,(0-^  (26) 

where  ^  =6.2-10  ^  U*  is  a  friction  velocity  . 

It  is  interesting  to  make  a  comparison  between  the  Toba  spectmm  and  the  Kolmogorov  one  (17)  or  (19).  In  order  to 
do  it  the  wave  energy  flux  (22)  should  be  defined.  It  includes  the  value  of  wave  energy  dissipation  in  a  high 
frequency  range,  which  is  unknown.  It  can  be  estimated,  supposing  that  the  wind  input  energy  flux  is  expended  to 
support  the  high  frequency  energy  flux  directed  to  a  high  frequency  range  of  wave  spectrum: 
n  W/- 


p  =  |  |/(co,(3)5((o,[3)t/(oJ[3 


(27) 


-31  a)„ 


The  wind  wave  input  energy  increment  /  can  be  estimated  using  a  traditional  approximation  [Komen  et  al.,  1994]  . 
It  is  defined  as  follows: 


/(C0,p)  =  max<^0;  0.252^(01 

1  Pw 


^  a  ^  ? 
- cos(p-p(/)-l 

c 


(28) 


where  Pa  and  Pw  is  a  density  of  air  and  water,  correspondingly. 

The  angular-frequency  spectrum  approximation  should  be  substituted  into  (22)  in  order  to  define  the  energy  flux 
(22).  Unfortunately  an  angular  distribution  of  the  Toba  spectrum  is  unknown.  Nevertheless  it  can  be  suggested  that  it 
is  rather  wide  in  high  frequency  range.  The  following  estimation  can  be  obtained  for  the  isotropic  frequency 
dependent  spectrum: 


P=  |/(a))5(co)£/co  =0.25-10-^  6gC/,  J 


CO' 


OD>g/28t/* 


2%Um 

g 


-1 


t/co 


(29) 


^ 28 f/* 


Thus,  the  energy  flux  can  be  estimated  as: 

P  =  0.125  10“^  5^  f/. 

Using  (19)  and  (25),  the  first  Kolmogorov  constant  can  be  estimated  as  follows: 


an=- 


5((0)  (if 
,4/3pl/3 


hgU* 


01g‘^l^9iWh  g 


=  2.175^/^  =0.33 


(30) 


(31) 


Comparing  this  value  with  that  obtained  with  the  help  of  d  irect  numerical  simulations  (see  Tabl.  1)  a  conclusion  can 
be  made  that  they  are  in  a  good  agreement.. 

Conclusions 


Direct  numerical  simulations  of  the  Hasselmann  kinetic  equation  for  gravity  waves  in  water  surface  confirms  basic 
predictions  of  the  weak-turbulent  theory.  The  kinetic  equation  for  surface  gravity  waves  is  investigated  numerically 
taking  into  account  an  external  generating  force  and  dissipation.  An  efficient  numerical  algorithm  for  simulating  non¬ 
linear  energy  transfer  is  used  to  solve  the  problem. 

Three  stages  of  wave  development  are  revealed:  unstable  wave  energy  growth  within  a  range  of  external  force 
impact,  fast  energy  spectmm  tail  formation  in  high  frequency  range  and  establishment  of  a  steady  state  spectrum.  In 
both  isotropic  and  non-isotropic  cases  the  spectra  are  found  out  to  be  dose  to  the  Zakharov-Filonenko  spectmm 

4 

CO  in  the  universal  range.  Reliable  estimations  of  the  Kolmogorov  constants  are  found  out.  Formation  of  this 
asymptotic  spectrum  happens  explosively.  Accurate  estimations  of  the  first  and  second  Kolmogorov  constants  are 
obtained. 

A  good  agreement  between  the  Toba  experimental  data  and  our  results  obtained  with  the  help  of  direct  numerical 
simulation  is  observed. 
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